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The author's gravitational potential for an oblate planet, expressible in terms of oblate spheroidal 
coordinates, can be generalized by means of a metric-preserving transformation of the associated 
Cartesian system. This preserves separability of the problem of orbital motion when the potential 
coefficients A. i and J■^ are taken into account. The inclusion of 7-2. i is not of practical importance, 
but has a clear |)liysical interpretation. The inclusion of 7.?. however, is of considerable practical 
importance, (>(^rmittin<; a more accurate treatment than that jiiven by perturbation theory. On the 
other hand, its physical significance is less clear, since including./.} depends on translating the origin 
of sph<*r<)idal c(K»rdinates by a distance ^= 2 J2^\J:\\ ecpiatorial radii. This distance, amounting to 
7 km in the case of the earth, is much greater than any change in the geoid produced by ./a. It is 
clearly related to the long-periodic terms of perturbation theory and turns out to be equal to the dis- 
placement by ./.{ of the plane of symmetry of those exactly elliptical polar orbits which are possible 
solutions with the spheroidal model. 
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1 . Introduction 

If one neglects the axial asymmetry of an oblate planet, the gravitational {)otential V outside 
it is given by 



r 






(1) 



Here /x is the product of the gravitational constant and the planet's mass, re is its equatorial radius, 
and r and 6 are the planetocentric distance and latitude of a field point. The PnS are the Legendre 
polynomials and the As are the coefficients of potential, with A ~ 10~^ and all higher J's ~ 10~^ 
or smaller for the earth. The first harmonic, /i= 1 , can be eliminated from (1) by taking the origin 
at the center of mass. 

If r is the position vector of a satellite of such a planet, its motion around the planet is described 
by r = — V F. In earher papers^ [V1959a, b] I discovered a choice oUn^ that reduces the solution 
of this differential equation to a separable problem. This choice is 

,12, = (- lf^\J^- J2k^i = (- iy\/i./^ (2) 
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corresponding to the spheroidal potential 

M(p + i7§) 

p^-fcV' ^^ 

p, 17, and if being spheroidal coordinates. The relations among rectangular, spherical, and 
spheroidal coordinates are 

x + ir= r cos ^ exp i(^ = [(p^-h c^)(l — r/^)]^/^ exp ic^ (4.1) 

z = r sin ^ = pT9. (4.2) 

Here 

c'^^rlJ. 3 = -reJi (5) 

and 8 is the distance, taken positive northward, from the origin to the center of mass. Since 
6 = when they coincide, this potential is symmetric with respect to the equatorial plane. Inver- 
sion of the separated equations led to an accurate reference orbit that takes into account the 
zeroth and second harmonic exactly and about two-thirds of the fourth harmonic [V1961a, b, 1962]. 
Shortly after the publication of V1959a, Brouwer and Pines [Brouwer and Clemence 1961, 
pp. 573, 592] found that the potential (3), with 8 = 0, could be discovered in another way. If one 
imagines half the mass of the planet placed on the polar axis at z = c' and half the mass at z = — c', 
the corresponding two-center potential is 

where 



■-x^^y^ + iz^c'f. (6.1) 



This corresponds to the separable problem of two fixed centers. If one now replaces c' by c V— 1 , 
the expression (6) remains real and leads to the separating potential (3) of the oblate planet. The 
significance of this procedure is traceable to the imaginary transformation that connects prolate 
and oblate spheroidal coordinates. 

Some time later Aksenov, Grebenikov, and Demin [1963] discovered that if /xu ^2. n, and r^ 
are all complex, with /Xir^' and yi-ir^^ conjugate, the potential 

F=-^-^ (7) 

also leads to separability and enables one to fit not only p. and Jz in (1) but also/;}, with the origin 
still located at the center of mass. An endeavor to understand this possibility in more physical 
terms led to the present paper. 

2. A Further Summary of Previous Results 

Let 

^^p/c, (8) 
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where c is the adjustable distance in (4.1). The potential (3) was the most ti;eneral function of 
oblate-spheroidal coordinates that satisfied certain conditions: 

a. It must be a solution of Laplace's equation. 

b. It must not be infinite everywhere on the polar axis. 

c. It must lead to separabihty of the Hamilton-Jacobi equation. 

Conditions (b) and (c) imply axial symmetry, so that no explicit assumption was really necessary 
for that. It happened, however, that a further condition was used, but not explicitly stated: 

d. The spheroidal coordinates are centered around the center of mass and around the polar 



3. Transformations of the Rectangular Coordinate System 

Giving up this fourth, unstated, condition permits an extension of the results. To see how, 
note that the Hamiltonian can be expressed as 

^^9 ;i7^ ^V(^.y^) (9) 



and Laplace's equation as 
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--(£-$-£)-»■ 



Postulating that V(^, j]) be a potential that satisfies (10) and leads to separabihty, we then ask what 
changes one can make in the rectangular coordinate system on which ^, 17, <p are based, without 
losing these properties. Since invariance of ds^ = dx' -\- dy^ -\- dz^ leads to invariance of V^, the 
answer comes easily: any metric-preserving transformation. Actually the properties are preserved 
even if ds^ and V'^ change by constant factors (conformal transformation). If, however, we write 
eqs (4) in the form 

x-hiy=c{{e^in''){l-r)'')Yl^' exp iip (11.1) 

z = c^V. (11.2) 

it is clear that we already have a disposable scale factor, viz, c, so that we cannot gain any added 
advantage by such a generahzation. Also we do not gain anything by considering reflections, 
since these do not introduce any extra parameters that might be used to fit more/'s. Hence our 
practical result: to obtain the greatest generalization of the results, we need consider only trans- 
lations and rotations of the rectangular coordinate system. 
Now from V 1959a we have the result that 

V=bo Re (f +ir?)-'+6, Im {^-hir))' (12) 

is the most general potential, based on the polar axis as 2-axis and on the center of mass as origin 
for the spheroidal coordinates, that satisfies conditions a, 6, and c. We may now translate and 
rotate the associated rectangular system and (12) will still have the desired properties a, 6, and c. 
Since a translation is described by three parameters and a rotation by three parameters, it would 
appear at first that we could fit six more potential coefficients. Actually the number is only five, 
since of the three rotation parameters, it takes only two to fix the new z-axis. Since the third is 
used to fix the rotation of the new x and y axes around the latter, it is immaterial with a potential 
like (12), which is symmetric around the new z-axis. 



4. Rotation 

We assume that the polar axis is not an axis of symmetry, but that there exists such an axis 
OZ in the planet, with angular coordinates i// and j8. A plane through O, perpendicular to OZ, 
intersects the plane Oxy in the line of nodes A^, which we take as the new Z-axis. Then if i, j, 
and k are unit vectors along Ox, Oy, and Oz, unit vectors k' and i' along OZ and OX are given by 

k' = i sin i// cos/3-(-j sin i// sin ^-h k cos i// (13.1) 

i' sin ip = klX k = i sin i//sin/3— j sin i//cos/3, (13.2) 



that 



i' = isinj8-jcos^ (13.3) 



j' = k' X i = i cos i// cos/3 + j cos i// sin^ — k sin i//. (13.4) 



Then since 



i:r+jy4-kz = i'A + jT+k'Z, (14) 



we obtain the equation 




(15.1) 



where A is the square rotation matrix 

fsin /3 cos ijj cos (3 sin i// cos ^\ 

— cos j6 cos i// sin ^ sini//sinj8 • (15.2) 

— sin i// cos i// / 

Here OZ is along the axis of symmetry of V= b(^e(^ -\- irj)~^ -\- bi Im(^-h ^17)"^ having angular co- 
ordinates i// and j8, the axis OX is along the line of nodes, and the sense of OY is given by the right- 
hand screw rule. 



5. Translation Plus Rotation 

If we also translate the origin from to O' by a vector 

00'^S = i5i +J52 + k53 (16) 

we find 

/x\ /si\ /X\ 

(17) 




111(1 



Y\ = A 



x — sA 


/sin^ 


y— 52 = 


= cosi//cosj8 


z — sj 


\sin \\} cos P 



-cos^ \/x-Sx\ 

C()si//sinj8 — sini|/ My — 52 
sini//sirij8 cos i// /\z — sj 



(17.1) 



Then 

Z= (jc — 5i) sin i// cos j8+ (y — 52) sin i// sin/3-h (z — 5.0 cos i//= (r — s) • k' (18) 

Also, if R is the position vector i'X-h j'F+k'Z, then R = r — s and 

R^-=-r^ + s'^-2^'v, (19) 

Our separating potential is then 

F=6oRe(^-hiT]) '-f/;, lm(^-h/r/)-', (20) 

where 

A' + iT=c[(f2_^l) (l-T/^)]>/-^expi(^', (21.1) 

Z-c^7/ (21.2) 

and (^' is the (^-angle relative to O'XYZ. 

6. Spherical Harmonic Expansion Relative to the Original System Oxyz 

In order to fit the coefficients of potential, we must now expand the potential V{^, rj) in spheri- 
cal harmonics relative to the original system Oxyz. Since we shall now obtain tesseral harmonics 
also, we have then to compare the expansion to be obtained with 



F=-^ 



1 — 2 ( — 1 JnPn (cos 0) + ^ ( — 1 2 P;;' (cos 6) (C„, ,/, COS nup + Sn, m sin mif) 



(22) 



In this section only, we use to mean colatitude, rather than latitude, of a field point relative to 
Oxyz. This is to facilitate use of the addition theorem for spherical harmonics and because i// is 
the colatitude of O'Z. 

To obtain the expansion of V{^, r\) in (20), we need to ex[)and (^-(-^17)"'. Now 



(^ + /T;)2=^2_^2_^2t^T; 



and from eqs (21) 



^2_^2^^_1 



R^ 

C2 



21^7] = 2i-' 
c 



(23) 

(24.1) 

(24.2) 



From (23), (24), (18), and (19) it then follows, if 8^ = 52, that 



(^+jr,)2 = - 



o2 _ p2 f. 2 

1+^— T^-2i-s-k'+4r-(ick'-s 
n r n 



(24) 



But 



(ich' — s)^ = s^ — 2ics -k' — c^, 



(25) 



since k'- = 1, so that 



i^-^ivr=-, 



(ick—sY 2 

1 H z ^ ~ Uck — s) • 1 ,■ 

r- r 



(26) 



where 1,- is a unit vector along r. Then 



(^+17})-' = - 



r \ r 



-1/2 



(27) 



Since c and | s | will both be small compared to r, we have 



-(ick — s)- 1, -f 

r \ r 



and we then have the valid binomial expansion 



<^ 1 



(28) 



i^^iv)-'- 



l-i(.ck'-s).l.-i(^^?^^ 
r 2 \ r 



(29) 



Now the term — cr -{ich'~s)'lr will give rise to first harmonic terms unless we choose s properly. 
From (20) we find that the choice 



Z>os — 6i ck' = 



(30) 



will ehminate terms in r~^ from F(^, rj) and thus place the origin at the center of mass. Thus 
we must choose s to be either parallel or antiparallel to k' (implying that the new Z-axis must pass 
through the old origin, even after the combined translation plus rotation), so that 



s=sk' (31.1), bos-bic = 0, 

where the sign of s is still undetermined. Then 

ick' - s = (ic - 5)k' (32.1), (ick' - s)^ = (ic - 5)- 
Insertion of (32) into (27) then gives 



{^+ivr'=- 



l--(5-ic)k'-l.+ f^— ^ 

r \ r 



-1/2 



Since 

Ir — i sin 6 cos ip + j sin sin c/? -f k cos 0^ 
it follows from (34) and (13.1) that, if U is the angle between k' and 1 r, 

k' • 1, = cos U= cos i// cos 0-\- sin i// sin 6 cos ((p ~ /S). 
6 



(31.2) 



(32.2) 



(33) 



(34) 



(35) 



Then, from (33) and (35), 



(^ + ^'^)-^=- 



1 (s — ic) cos U-\- 



r 

s — ic 



s — ic 



-1/2 



;2 V^M'»(-f/). 



(36) 

(37) 



Insertion of the addition theorem for spherical harmonics 

^ {n — m) 

■'j {n -h m) 



Pnicos U) = Pnicos ijj) Pn(cos ^) + 2 £ j PX'(cos i//) P;f(cos 6) cos(m(^ - 7^/3) (38) 



into (37) then yields 

/s — ic 



i^+iv) 






PnicOS l//) P.(C0S 0) 

+ 2 y '—^ 4 P;;'(c()s ijj) P;;'(cos 6/) cos(m<^-m/3) 

.r', (^+^)! 



(39) 



With apphcation of (20) and (31.2) to (39) we may readily verify that the terms involving 7], Ci, i, 
and Si, 1 are missing from (22), so that we have indeed placed the origin at the center of mass. 
We began with eight disposable parameters bo, b\, c, s\, s^, 53, i//, and ^ and after ehminating the 
first harmonics we are left with five, viz, 60, 5, c, i//, and /3, b\ being fixed by (31.2). 

On comparing (39) and (20) with (22), we next find that to fit the zeroth harmonic, we need 



Then 



Let us now put 



b()C = — fji. 



b() = and bi= r 



S = OTr cos 7 C = OTe sm y. 



so that 

s — ic = are exp (— iy)^ {s — ic)" = cr"r'/. exp (— iny). 
From (39), (20), and (22), we then find for the zonal harmonics, 

fJir'/Jn = cPnicos i//) [bo Re(5 - ic) " + 61 lm(5 - ic) "] . 

Application of (43) to (44) then gives 

sin in— l)y 



sm y 



Jn = O-"Pn(c0S ijj) 

yielding the expected results 70 = — 1 , Ji = 0, and 

J 2 =- 0-- P2(cos i//) 

J'.i = 2a'-^ P'.](cos i//) cos y 

7 



(40) 
(41) 

(42) 
(43) 

(44) 
(45) 

(46.1) 
(46.2) 



, . r^ , , X sin 3 y 

J^ = 0-4 P4 cos i^ —. ^ • (46.3) 

sin y ^ ^ 



For the tesseral and sectorial harmonics we find 

Cn,m^iSn, m = -2 J^^^^ cx" P-(cos i//) ''" ^^^ ~ ^^^ exp(im ^8), 
{n-\-m)l sin 7 

leading to the expected results Ci, i = Si, i = and to 



(47) 



C2, 1 = - y sin 2i// cos ^ (48.1) 

0-2 
S2, 1 = - y sin 2i// sin /3 (48.2) 

C2,2 = --rsin2i// cos2)8 (48.3) 

S2,2==-xsin''i//sin2j8. (48.4) 

We now have four parameters remaining, cr, i//, 7, and j8, to be fitted to certain of the potential 
coefficients above. Clearly we must fit to J2- The remaining choices are to fit to J3 or J 4 and to 
Cz, 1 and S2, 1 or to Cz, 2 and Sz, 2. Now it is clear that in any practical application i// would have to 
be a very small angle. Let 

J2,i = {ci,^si,yi^ 72,2= (Ci, 2+ Si, 2)^/2. (49) 

From (48), Jz, z would be of order ./I, i, making Jz, i much larger than Jz, z- This is contrary to all 
experience, since for the earth ^2,2 is a measure of equatorial ellipticity and72, i of the variation 
in latitude produced by wobbling of the polar axis. We should therefore fit not to C2, 2 and S2, 2 
but to Cz, 1 and S2, 1. The resulting much too small values for Cz, 2 and ^2,2 simply mean that we 
cannot fit them by this model. 

7. The Physical Significance of the Rotation ijj, (3 

As we have defined it, ip is the angle between the polar axis and the axis of symmetry of the 
separating potential V=bo Re{^ -\- ir})'^ -\- bi lm{^ -\- irj)~^ . Now nonvanishing of i// leads to non- 
vanishing of ^2, 1, which would actually vanish if Oz were a principal axis of the primary body. It 
is thus natural to identify the axis of symmetry of V{^, 7]) with that principal axis which is closest 
to the axis of rotation. 

For the earth i// ~ 0.1", the mean variation of latitude associated with the wandering of the 
pole. For the moon it is likewise expected to be very small. 

Furthermore, when we use a potential that has symmetry about some axis other than that of 
rotation, the potential, and thus the Hamiltonian, become explicit functions of time, if we use an 
inertial reference system. This circumstance would defeat the separability. There is another 
possibility, that of letting the system Oxyz be rigidly attached to the primary body, as we have actu- 
ally done above. We should then have to treat the effect of rotation on the orbit of a satellite by 
means of apparent forces. They are not explicit functions of time, but they also destroy the 
separability, unless they are so small that we could neglect them. There might be some point in 
such a procedure for the moon, since its rotation is relatively slow, but even for the moon the 
apparent forces will undoubtedly be a good deal larger than those arising from Jz, i- If we could 
have fitted 72,2, the story might have been different. 



All in all the rotational invariance of F(^, 17) does not appear to lead to any practical results. 
1 have included a treatment of it partly to verify that we do not miss anything useful hy abstaining 
from its use, but mostly because its physical significance is so transparently clear. When we make 
such a rotational transformation, the new axis that we find is an axis of symmetry of the resulting 
potential. This stimulates us to search for a physical meaning of the translation s, wfiicli turns 
out to be of real practical value, but whose physical significance is much less apparent. 



8. Results of the Translatlonal Invariance of V{^, rj) 

On putting i// = 0, we thereby ehminate the rotational transformation and obtain the results 
of the pure translational invariance. We then have left only two paramters, a and 7, to be fitted 



to the Jn's. By (45) we thus obtain 



_ 8in(//-l)y 
sm y 



(50) 



leading to 

./2-0-- (51.1) 

J:i = 2(T'' cosy (51.2) 

./4 = a-4(4 cos^y-1) (51.3) 

J^ = 4'(T'' cos 7(2 C0S-7-I). (51.4) 

We must certainly fit to ./j- Thus 

cT=Jll^' (52) 

Since (51.3) and (51.5) lead to 

4 COS- 7=^ — 7^ — , (52.1) 

J2 

we could fit to J4 only if ,/4+,/| were positive. Although the various investigators still find widely 
differing values of ^4, they now agree that J^^Jl is negative. Thus we cannot fit to J4. 
On fitting to ./.{, we then find 



cosy=i-JJ2-'l'\ (53) 



which leads by (42) and (52) t( 



^=^/. j^ (54.1) c = rj^^ {^~IW ^^"^'^^ 

With ./2 = (1.083)10-3 and A - - (2.4)10-«, we find j JU^'' - 4 X 10-^, so that the new value for c 

4 

is about 2 parts in a thousand smaller than the value arrived at in V1959a. With i|/ = 0, we have 

s = 5k, where 5 ~ — 7.09 km. The values of J4 and J^ corresponding to the above values of cr and 

7 are 



-Ji + jJ (55) 



J5 = -2J:Jj2--jj- (56) 

The value of /4 given by the model is thus little changed from the earlier value of — y| (V1959b) 
and the values of J5 and all higher J's are negligibly small. The fit to the potential is thus exact 
through the third zonal harmonic and accounts for at least two-thirds of the fourth zonal harmonic. 

9. The New Form for the Potential 

From (20) and (41) and p = c^, we now have 






y—'Tr:^ (57) 



If we place 

5 = -S, (58) 

where 

8=5/v^ (59) 

is the magnitude of the distance from the new origin of spheroidal coordinates to the origin, we 
obtain 

y=-t^-4~fi- (60) 

This of the same form as eq (49) in V1959b and 8 has the same meaning that it had there, viz, 
the distance northward from the origin of spheroidal coordinates to the center of mass. The differ- 
ence is that in the earHer paper the origins of spheroidal coordinates and of spherical harmonics 
were always the same. Had I kept 8 7^ but expanded the spheroidal potential in spherical 
harmonics about the center of mass, I would have been able to fit ^3 at that time. However, the 
more general expressions for the quartic G{r\) that appear in eq (63) of the earher paper can now 
be used. 

The connections between rectangular, spherical, and spheroidal coordinates are now of course 
a little different. They become 

x^iy—r cos 6 exp i(p= [(p^ + c^) (1— 17^)]^^^ exp i(p (61.1) 

z = r sine = -di-pr] (61.2) 

Actually the only difference is the appearance of 8 in (61.2). 

Previously, with 8 = 0, a plane through the origin of spheroidal coordinates, perpendicular 
to the polar axis, had two important properties. It passed through the center of mass and was a 
plane of symmetry of the potential. Now, with 8 7^ 0, by (60) and (61.2), such a plane has neither of 
these properties. In contrast with the rotational transformation, which was to an axis of sym- 
metry of the potential, the translation s = — k8, although it has the merit of enabling us to fit Js 
also, seems difficult to understand physically. We shall find later that a solution for some of the 
possible orbits will give us some physical grasp of the translation. 
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10. Remarks About Resonance Denominators 

Before solving for the orbits, it is desirable to discuss briefly the bearing that the present topic 
has on the singularities found in perturbation theories of satellite orbits. If / is the inclination, 
the resonance denominator 1-5 cos^ / has been found, by Brouwer 1959, Garfinkel 1959, Kozai 
1962, V1963, and others, to occur with numerators 74+71 and J^, but not with the numerator J3. 
In other words they are absent li J^^J\ and the odd /'s vanish. One could previously understand 
their absence in such a case, because it corresponds, by (2), to the author's earlier spheroidal 
potential (3), with 6 = 0. Such a potential, leading to separability, is not expected to lead to such 
resonance denominators in a perturbation theory. E.g., the solution of V1961b, V1962, gives exact 
secular terms and periodic terms which can be worked out to any degree of accuracy, with every- 
thing remaining finite when 1-5 cos^ 7 = 0. Thus Kozai 1962 was able to do a partial check of his 
results by seeing whether his resonance denominators all vanished ior J.i = J^=. . .=0 and 
A + 71 = 0. 

We now have a more stringent criterion for checking the accuracy of long-periodic terms in 
perturbation theories. The resonance denominators should all disappear when the ./'s are given 
by (50), with a = ,/-/' and cos y= 2J'aJ2^''-. I have tried to check Kozai [1962| in this way, but find 
the going difficrult. Despite the great complexity of that paper, it may be that he has not carried 
terms to a sufficiently high order to make the test applicable. 



1 1 . Remarks About the Translation Parameter 8 

Since 8= 2 nJo^lJ-.]] ~ 7 km for the earth, it is indeed a large quantity to be j)rodu('ed by the 
nonvanishing of J^. The latter produces a maximum change in the geoid of only about 20 m, a 
quantity of the order re\J;i\. To seek for some physical significance of 6, I shall consider the 
effect that Js has on some particularly simple orbits, viz, orbits which would reduce to circular 
equatorial orbits and circular polar orbits for V= — fxr' K 

1 2. The Special Orbits 

To find these special orbits, we can take over some results from V1959b and V1961a, without 
re-solving the whole problem. The generahzed momenta pp and pv are given by 

(62.1) 

(62.2) 

where F{p) is a quartic with physical zeros pz^ p\> and G{y]) is the quartic 

G[r]) = - dr^^ {\ — 7)-) [di^ 2^7)8 + 2axc~y)'-) (63.1) 

= di — al-\-2ij.r]8 + (2a, r- — ai)7]- — 2/>tT7"^8 — 2aic-r]^- (63.2) 

In eqs (63) a, is the energy, a-^ the z-component of angular momentum, and az the separation 
parameter, reducing to the total angular momentum in the Keplerian case. The physical zeros 
of G{r]) are r]o and 171, with 

1^170^171^-1, (64) 

since r}'^ ^ 1. 
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' p^ + r^ 


f^- 


p'--^ r- 




V^ 





We introduce the definitions 



ao= — -: — , eo = IH ^^ , io= cos ^ (a3/<^2), Po = ao(l — e^) 

a - *(pi + p2), e ^ ^^^^, p - a(l -6^). (65) 

P2-^Pl 

Then for an orbit with e = 0, F{p) remains zero, and p is constant, by (62.1): 

p = a = p. (66) 

From (4.14) and (4.15) of V1961a, we find for an equatorial orbit that 

aopo=d' -h 3c2 + OUl) (67.1) 

ao=all^^Wo{J'i). (67.2) 

For a polar orbit, it is clear from (61.1) that 17^ reaches the value unity, so that 

T7o = +l i7i=-l- (68) 

From the following paper we also have for a polar orbit 

e = €o a = ao p = Po. (69) 

If e = 0, then eo = 0, so that ao = po and 

p = a= ao = po = p. (70) 

13. An "Equatoriar Orbit 

For the case 7.3 = and thus 6=-0, it follows from V1961a that there is a solution corresponding 
to rjo= 7)1 = 0, so that r) then remains zero and the orbit is confined to the equatorial plane 2 = 0. 
For J3 ^ and thus 6 7^ 0, such an orbit does not exist: a satellite moving in the plane z = would 
be driven out of that plane by the force 



Fs = V 



pi'^Psi^ine) 



r 



(71) 



arising from the third harmonic. 

When 8 7^ 0, we see from (61) that the only circular orbit that always lies in a plane perpendicu- 
lar to the polar axis is given by 

p= constant 17 = constant. (72) 

But p can be constant only if pi = p2 and then 

P = «- (73) 

The coordinate tj can be constant only if 

Vo = Vi = v'- (74) 
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To find the z-coordinate of this plane we must therefore find iq' . Now when 170 = 171 = 17', 17' is a 
double zero of G(r)) and is therefore a solution of 

Thus, by (63.2) and (75), 17' is a solution of 

2/x8 + 2{2ai(^ - ai)r]' - 6fjLr}"'d - Saic'r]"' = 0. (76) 

if 8=0, we know that 17' = 0, so that we may now assume 17' to be small. On using the notation 
in (65) and solving (76) by successive approximations, we then find 



V 



36^ 

Pi. 



+ 0{Jt)' (77) 



Then, by (61.2), (73), and (77), the 2-coordinate of the plane of this "equatorial" orbit is 

«8 



po+c~/ao 
From (67.1) and (67.2), we now find 



+ OUt)' (78) 



that 



Po+ — = a+'-^+0{fi), (79) 



2'=-a+ , , .f .,. ., + (Kn) (80) 



z' = -Y^^,^ + (Xm (81) 

I + 6c-ia- 

= -^6 + W;S), (82) 

a- ~ 

since 8 = 2r^j7'|./:}| is of order Jz- By (54.2) we now have c^ = r^ 2 ~^ (Xfi) , so tliat finally 

2' = |^ya + 0(/i)- (83) 

For close orbits this gives z' ~ — 22 m, a reasonable figure that is of the same order as the maximum 
change in the geoid produced by J^. 

As a check on the above expression, one may apply the perturbation theory for Jn given in 
V1963, as applied to an orbit for which e = and sin / = 0. Since these quantities appear in de- 
nominators of the expressions for the variations of the Delaunay variables, one has to carry out 
limiting processes. It turns out that the long-periodic variations give 8z = 0, but that the short 

3 
periodic ones give 8z= - r|a~V3, equal to the above z' . Thus the present result is in agreement 

with perturbation calculations. 
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The result can also be verified directly from the equation z = — dVldz. Here 

-F=^ + iMr|i.(i-^)-fi^r^3(|-f) + . . .. (84) 



The 



z=- — + 2 /x^^2 (- 7 + -TTJ + ^f^^V'^ [;7.--;r^ ^) + • • •• (^5) 

For an orbit of constant z and constant radius a, we have 2 = and 

r^ = a' + z' = a' (l + ^)- (86) 



On neglecting powers of z higher than the first, we then find 

_ 3 j^Js _ 3 r? 



a- + - ry 2 



Js^OUl). (87) 



as we found before. 



14. A Polar Orbit 

Since we cannot interpret the translation parameter 8 as the displacement produced by Js 
of the simplest equatorial orbit, let us now try the simplest polar orbit. 
We first consider the case 6 = 0. Then by eqs (61) 

;c- + r = (p- + c^)(l-T7-) z = piq. (88) 

If e = 0, then p= a and 

leading to 

4^ + 4 = 1- (89) 

This is the equation of the coordinate spheroid p = a and simply means that when e = 0, the satel- 
hte's motion is restricted to it. If the orbit is polar, it is then the intersection of a meridian plane 
with such a spheroid. It is thus an ellipse with semiminor axis a, semimajor axis (a^ + c^)^/^, 
with center (not focus!) at the center of mass of the planet. It is not a Keplerian ellipse, since the 
law of areas does not hold; if R is the distance from a focus to the satellite, it could not be produced 
by a potential — piR'^ acting on the latter. That my original spheroidal potential [V1959a] leads to 
the existence of such an orbit was first noticed by Aksenov, Grebenikov, and Demin [1964]. 
By (88), this elliptic polar orbit lies between the limits 

z = ±a. (90) 
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If we now consider the case 8 7^ and still choose e = 0, then p = a as before and by eqs (61) the 
orl)it is the intersection of a meridian plane with the spheroid 

-^T^^- 7^ = 1. (91) 

cr -T c- a- 

It is thus an ellipse unchanged in size or shape, but with z now varying between the limits 

Z2 = — ^ -\- a 

zi = -8-a, (92) 

by (61.2). 

For J3 = 0, the plane of symmetry of all such elliptic polar orbits was the plane z = i). For 
J:i # 0, their plane of symmetry is now 

2 = -8= 2-/v./,7\/.v (93) 

By (65) and (70), the energy ai=— 2 ixjao and ^0 = a, so that all such orbits of f\xed a have the 
same energy. From the conservation of energy one can show that any such orbit can be produced 
by first placing the sateUite in the plane of symmetry, at a distance ro from the polar axis, and 
tiien projecting it parallel to the polar axis with initial velocity 

Zo=(fJilayl'\ (94) 

where 

<r = ii-c'- (95) 



It appears that we have now 


arrived at a sort of physical interpretation of the translation 


param- 


eter 8. When we go from the separable mod 


el with J:i = 


to the improved one that fits 


J,u the 


necessary translation of the origin of spheroidal coordinates 


is equal to the displacement 


of the 


plane of symmetry of the elliptic 


polar orbits. 








A certain verihcation can be obtained from the perturbat 


ion theory of V1963, by placir 


ge = 


and 170— 1 into tlie long-periodic 


terms there gi 


ven. We hnd 


for the variation of z 




dz = 


1 J.i . ., , 

2 -r re sm- i//- 
J2 


J, ... 

-y /-p sin / sin 


^ cos l//. 


(96) 


Thus at either pole, with sin i// = 


±1, 










3z = 


1 .h 




(97) 


corresponding exactly to the above displacement. Moreover 


short-periodic averaging gives 






+ A> sin^ ^ 




(98) 


and long-periodic averaging then 


gives 










8~z = 


1 7,3 




(99) 


the same displacement as above 












\ 
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We may sum up the physical results as follows. If we consider the equation r = — V V, 
where V includes only even zonal harmonics and where all Jk's higher than Jo are much smaller, 
there must exist a class of polar orbits which are close to true ellipses, viz, those polar ellipses 
corresponding to A — — J|, Jq^JI, etc. All these true ellipses of a given energy generate a sphe- 
roid, and varying the energy furnishes the whole set of coordinate spheroids used in the sphe- 
roidal method. The equatorial plane is the plane of symmetry of these spheroids. It is also a 
plane through the center of mass and the plane of symmetry of the potential when only even zonal 
harmonics are included. 

When we include J^ also, if all higher zonal harmonics are small, there must again exist polar 
orbits which are very close to true ellipses, viz, those polar ehipses corresponding to the special 
choices made for the higher harmonics in the present paper. Again these new ellipses generate 
the new coordinate spheroids. The plane through z = — 6, perpendicular to the polar axis, now 
does not pass through the center of mass and is now not a plane of symmetry of the potential, 
but it is again the plane of symmetry of the elliptical polar orbits. 
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